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Abstract-Based on a new version of the Hipparcos catalog and currently available radial 
velocity data, we have searched for stars that either have encountered or will encounter the 
solar neighborhood within less than 3 pc in the time interval from —2 Myr to +2 Myr. Nine 
new candidates within 30 pc of the Sun have been found. To construct the stellar orbits 
relative to the solar orbit, we have used the epicyclic approximation. We show that, given 
the errors in the observational data, the probability that the well-known star HIP 89 825 
(GL 710) encountering with the Sun most closely falls into the Oort cloud is 0.86 in the time 
interval 1.45 ± 0.06 Myr. This star also has a nonzero probability, 1 x 10~^, of falling into 
the region d < 1000 AU, where its influence on Kuiper Belt objects becomes possible. 

DOI:10.1134/S1063773710030060 

INTRODUCTION 

Interest in the problem of close encounters of field stars with the Solar system stems from 
the fact that the passage of a star can lead to various kinds of perturbations of Solar system 
objects. 

Thus, for example, the interaction of a star with the Oort comet cloud can give rise to 
comet showers reaching the region of the major planets (Hills 1981; Weissman 1996). The 
formation time scale of a comet shower is ~1 Myr after the passage of a star. According 
to Weissman (1996), the Oort cloud is a spheroid with a semimajor axis of about 10^ AU 
(^0.5 pc, since 1 pc=206 265 AU) elongated toward the Galactic center and a semiminor 
axis of about 8 x 10^ AU. The question about close encounters of stars with the Sun within 
r < 2 — 5 pc was considered by Revina (1988), Matthews (1994), and MuUari and Orlov 
(1996) using various ground-based observations and by Garcia-Sanchez et al. (1999, 2001) 
based on Hipparcos (1997) data in combination with stellar radial velocity data. As a result, 
156 Hipparcos stars from the solar neighborhood 50 pc in radius that either have encountered 
or will encounter with the Solar system within r < 5 pc in the time interval ±10 Myr are 
known to date (Garcia-Sanchez et al. 2001). Having analyzed these data, Garcia-Sanchez 
et al. (2001) estimated the frequency of close (within 1 pc) encounters of stars with the Sun 
to be 2.3 ± 0.2 encounters per Myr and, after a correction for the Hipparcos incompleteness, 
this value increases to 11.7 ± 1.3 encounters per Myr. 

Analyzing the possibility of even closer encounters is of current interest. For example, 
Kenyon and Bromley (2004) argue that only the passage of a star can explain the peculiarities 
of the orbit of the minor planet 2003 VB12 (Sedna). The various kinds of influences of such 
a star on Kuiper Belt objects were shown to manifest themselves at characteristic distances 



r < 1000 AU (~0.005 pc). The Kuiper Belt proper extends from the Sun to a boundary of 
Ri50 AU. 

The goal of this study is to search for candidate stars closely encountering with the 
Sun based on a new version of the Hipparcos catalog (van Lccuwcn, 2007) and currently 
available stellar radial velocity data. We solve the problem of statistical simulations by 
taking into account the random errors in the input data and estimate the probability of a 
star penetrating into the Oort cloud region and into the region of a possible influence of the 
star on Kuiper Belt objects. 



1 THE DATA 

We use stars from the Hipparcos catalog (ESA 1997) while taking the new proper motions 
and parallaxes from a revised version of this catalog (van Leeuwen, 2007); the stellar radial 
velocities are taken from the Pulkovo Compilation of Radial Velocities (PCRV) (Gontcharov 
2006) created as a result of implementing the OSACA project (Bobylev et al., 2006) and 
containing radial velocity data for about 35 000 stars. Note that, in contrast to the larger 
CRVAD-2 catalog of radial velocities (Kharchenko et al., 2007), the PCRV catalog contains 
only the stars with random errors in their radial velocities within 10 km s~^. 



2 THE METHODS 

2.1 Orbit Construction 

Wc use a rectangular Galactic coordinate system with the axes directed away from the 
observer toward the Galactic center (/ = 0°, 6 = 0°, the X axis), in the direction of Galactic 
rotation {I = 90°, 6 = 0°, the Y axis), and toward the North Pole (6 = 90°, the Z axis). 
The corresponding space velocity components of an object U, V, W are also directed along 
the X, y, Z axes. The epicyclic approximation (Lindblad, 1927) allows the stellar orbits to 
be constructed in a coordinate system rotating around the Galactic center. We apply this 
method in the form given in Fuchs et al. (2006): 

X{t) - X(0) + ^ sin(«i)+ 

+^(l-cos(«:t)), (1) 

F(t) = F(0) + 2yl(x(0) + ^)t- 

-^VOS) siniKt) + ^C/(0)(1 - cosUt)), 

Z{t) = sm{vt) + Z(0) cos{vt), 

where t is the time in Myr (we proceed from the ratio pc/Myr — 0.978 km s~^); A and B 
are the Oort constants; k — ^^/—A.Q,QB is the epicyclic frequency; is the angular velocity 
of Galactic rotation for the local standard of rest, VLq = A — B-^ v = \/ AiiGpo is the vertical 
oscillation frequency, where G is the gravitational constant and po is the star density in the 
solar neighborhood. 



The parameters X(0), F(0), Z(0) and U{0),V{0),W{0) in the system of equations (1) 
denote the current stellar positions and velocities. The displacement of the Sun from the 
Galactic plane is taken to be Z{0) — 17 pc (Joshi 2007). We calculate the velocities U, V, W 
relative to the local standard of rest (LSR) with {U,V,W)lsr = (10.00,5.25,7.17) km s"^ 
(Dehnen and Binney 1998). 

We took po = O.IMq/pc^ (Holmberg and Flinn 2004), which gives z/ = 74 km s~^ kpc~^. 
We used the following Oort constants found by Bobylev et al. (2008): A — 15.5 ± 0.3 km 

kpc~^ and B — —12.2 ± 0.7 km s~^ kpc~^; k, — 37 km s~^ kpc~^ corresponds to these 
values. 

Note that we neglect the gravitational interaction between the star and the Sun. 

As was shown by Mullari and Orlov (1996), even the simplest linear approximation 
provides a sufficient accuracy in the time interval from —2 Myr to +2 Myr compared to 
the method of constructing the relative orbits (the star relative to the Sun) in the Galactic 
potential. We choose the epicyclic approximation, because wc arc planning to perform 
statistical simulations with the computation of hundreds of thousands of orbits, while the 
application of more complex methods requires a huge computational time. 

2.2 Statistical Simulations 

In accordance with the method of Monte Carlo statistical simulations, we compute a set of 
orbits for each object by taking into account the random errors in the input data. For each 
star, we compute the encounter parameter, the shortest distance between the stellar and 



solar orbits d^nin — jAX^{t) + AY^{t) + AZ^{t). The stellar parameters are assumed to be 



distributed normally with a dispersion a. We add the errors to the equatorial coordinates, 
proper motion components, parallax, and radial velocity of the star using the nominal errors. 



First, we considered the solar neighborhood 30 pc in radius using ^35 000 stars with known 
space velocities. For each star, we constructed the orbit relative to the Sun in the time 
interval from —2 Myr to +2 Myr. In addition to the stars from the list by Garcia-Sanchez 
et al. (2001), we found several other Hipparcos stars. Data on the nine stars found with an 
encounter parameter c? < 3 pc are presented in Table 1; Fig. la shows the trajectories for 
six of these stars. 

Figure lb shows the model trajectories of four stars computed by taking into account 
the random errors in the observational data; 300 realizations were obtained for each of the 
three stars indicated in the figure. 

Table 2 gives characteristics of the stars common to the list by Garcia-Sanchez et al. 
(2001) that have the closest {d < 2 pc) encounters with the Sun. 

Table 3 gives characteristics of the common hst of stars that have the closest {d < 
2 pc) encounters with the Sun. Note that this is currently the most complete list with 
the restrictions that we used. In contrast to Tables 1 and 2, Table 3 gives the encounter 
parameters dmin and tmin with their errors estimated through statistical simulations. 

Table 3 also contains two stars from Table 1. The first of them is the well-known white 
dwarf WD 0310.688 (HIP 14754). Note that the error in its radial velocity (Pauh et al. 
2006) is largest among all of the stars listed in Tables 1 and 2. The data for the second 




3 RESULTS 



star, HIP 27887, arc very reliable. It is contained in the catalog of F and G dwarfs of the 
Geneva-Copenhagen survey (Nordstrom et al. 2004; Holmberg et al. 2007). 

Note that the compilations by Gontcharov (2006) and Kharchenko et al. (2007) give an 
obsolete radial velocity for HIP 3829, Vr = 263 ± 5 km s~^, taken from the GCRV (Wilson 
1953). A very close encounter with the Sun is obtained with this value: dmin = 0.96 pc at 
tmin = —16 thousand years. However, this star is the well-known white dwarf WD 0046-1-051 
(van Maanens star, vMa 2). According to more recent observations, its radial velocity is 
Vr = 54 km s~^ (Greenstein and Trimble 1967) or Vr{LSR) = -41 km s"^ (Aannestad 
et al. 1993). In this case, the encounter is considerably less close, d^nin > 3 pc (see also 
Garcia-Sanchez et al. 1999, 2001). 

Figure 2 shows the trajectories of HIP 89825 (GL 710) relative to the Sun computed by 
taking into account the random errors in the observational data. 

We made only 300 realizations for a clearer graphical presentation. It is obvious from 
the figure that very close encounters with the Solar system are possible for this star. More 
detailed simulations show that, for example, out of the one million model orbits found for 
GL 710, dmin < 0.48 pc in 855 902 cases and, hence, the probability of its falling into the Oort 
cloud is Pi = 0.86 and dmin < 0.005 pc in 94 cases and, hence, the probability of the star 
approaching the boundary of its possible influence on Kuiper Belt objects is P2 — 0.94 x 10~^. 

Figure 3 show the distribution of encounter parameters of HIP 89825 (GL 710) with the 
Sun obtained using one million model orbits. 

4 DISCUSSION 

(1) Among the Hipparcos stars we consider, there are nine stars common to the list by 
MuUari and Orlov (1996). The encounter parameters for some of the stars found by MuUari 
and Orlov (1996) are: dmin — 1-35 pc, tmin = 1050 thousand years for GL 710; dmin = 1-64 
pc, tmin = —529 thousand years for GL 208; and dmin = 1-89 pc, tmin = 88-7 thousand 
years for GL 860 A. The encounter parameters agree well with our results (Table2). There 
is only one exception, GL 710, for which the new data give a closer encounter, which was 
also pointed out by Garcia-Sanchez et al. (1999). 

Note that the star HIP 99461 (GL 783) we selected is present in the list by MuUari and 
Orlov (1996) but is absent in the lists by Garcia-Sanchez et al. (1999, 2001). 

(2) For an overwhelming majority of the stars, the encounter parameters derived here 
agree well with those from Garcia-Sanchez et al. (2001). 

The differences are most noticeable for the star GL 217.1, more specifically, |A(ij„j„| = 
0.38 pc and |Atmm| = 182 thousand years. According to Garcia-Sanchez et al. (2001), 
who adopted the radial velocity Vr = 20.0 ± 3.7 km s~^ that differs markedly from our 
value (Table 2), the encounter parameters for GL 217.1 are: dmin = 1-65 ± 0.27 pc and 
tmin = —(1046 it 163) thousand years. 

Computing the model epicychc orbits for the star GL 217.1 using the initial data from 
Garcia-Sanchez et al. (2001) with the adopted radial velocity error — 2 km s"^ yields 
dmin = 1.61 ±0.20 pc and tmin = — (1066±116) thousand years, and these parameters almost 
coincide with those from Garcia-Sanchez et al. (2001). 

This leads us to conclude that the existing discrepancy is undoubtedly related to the 
difference in input data. 

As we see from Table 3, GL 217.1 is the most massive, M — 2Mq, star among our stars. 



Its improved encounter parameters, dmin = 1-28 =b 0.06 pc and tmin = —(861 ± 40) thousand 
years, make it attractive for studying the close passages that could take place in the past. 

(3) The star GL 710 is of greatest interest to us, because it has a record close encounter 
with the Sun among all of the candidates known to date. According to Garcia-Sanchcz ct al. 
(2001), the encounter parameters of GL 710 are: dmin = 0.34±0.18 pc and tmin = (1358±41) 
thousand years, in good agreement with those we found (Table 3). 

GL 710 is the only star with a significant difference between the data of Tables 2 and 
3. This suggests that the conditions for its encounter with the Sun make the stellar orbit 
very sensitive to small changes in such initial data as the parallax and radial velocity. We 
clearly see from Fig. 3a that the distribution of minimum encounter distance dmin differs 
from a Gaussian one, the mode of this distribution is dmin = 0.27 pc, which is lower than 
the median value of dmin given in Table 3, while the distribution of encounter times is nearly 
Gaussian (Fig. 3b). 

(4) Our simulations show that among the candidates considered, only the star GL 710 
has a high probability of penetrating into the Oort cloud region. 

5 CONCLUSIONS 

Based on currently available space velocity data for about 35 000 Hipparcos stars, we searched 
for stars closely encountering with the Solar system. For this purpose, we took stars within 
30 pc of the Sun and determined their orbits relative to the Sun based of the epicyclic 
approximation in the time interval from —2 Myr to +2 Myr. We found nine new candidates 
in addition to the well-known list of such stars (Garcia-Sanchez et al. 2001). The use of the 
PCRV (Gontcharov 2006) and improved Hipparcos stellar parallaxes (Leeuwen 2007) made 
this possible. 

As a result, wc compiled the currently most complete (given the restrictions) list of 
Hipparcos stars that have close {d < 2 pc) encounters with the Sun. 

For the star GL 217.1, a well-known candidate for a passage close to the Sun, the new 
observational data were shown to change noticeably its previously known encounter param- 
eters with the Sun. The encounter parameters found here are: dmin — 1-28 ± 0.06 pc and 
tmin = —(861 ± 40) thousand years. Improving the radial velocity for the white dwarf WD 
0310-688 (HIP 14754) whose orbit passed at a distance dmin = 1-61 ±0.19 pc from the solar 
orbit about 300 thousand years ago is of current interest. Our statistical simulations showed 
that the star GL 710 has not only a high probability of penetrating into the Oort cloud. 
Pi — 0.86, but also a nonzero probabihty, P2 — 1 x 10~^, of penetrating into the region 
where the influence of the passing star on Kuiper Belt objects is signiflcant. 
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Table 1: Hipparcos stars within 30 pc of the Sun found here 



HIP 


SP 


dmin-} 

pc 


10^ yr 


7r± e^, 
mas 


K ± evr, 
km s~^ 


n 


km s^^ 


1 


2 


3 


4 


5 


6 


7 


8 


14754 


DA 


1.60 


-296 


97.66 ± 1.85 


33.8 ±3.2 


1 




27887 


K3V 


1.97 


-413 


76.95 ± 0.37 


30.7 ±0.3 


1 




99461* 


K2V 


2.05 


41.1 


166.25 ± 0.27 


-129.0 ±0.2 


2 


0.8 


24186 


MOV 


2.15 


-10.8 


255.66 ±0.91 


244.4 ± 0.2 


3 


1.1 


105766 


G5 


2.20 


341 


37.87 ±0.42 


-76.6 ±0.3 


2 


0.2 


30344 


KOV 


2.30 


-2035 


34.10 ±0.62 


14.6 ±0.3 


1 




26373 


KOV 


2.33 


-762 


39.82 ± 1.36 


32.5 ±0.3 


1 




104214 


K5V 


2.79 


18.6 


286.82 ± 6.78 


-65.9 ±0.1 


5 


0.7 


101217 


K7V 


2.80 


19.6 


285.88 ± 0.54 


-61.3 ± 0.1 


3 


0.3 



Note, n is the number of catalogs used to calculate the mean Vr] £Vr is the radial velocity error 
calculated from the external convergence (Gontcharov 2006); (*) HIP 99461 (GL 783) is present 
in the list by MuUari and Orlov (1996) but is absent in the lists by Garcia-Sanchez et al. (1999, 
2001). 



Table 2: Data on the known Hipparcos stars from the solar neighborhood 30 pc in radius encoun- 
tering with the solar orbit within d < 2 pc 



HIP 








7r± e^r, 


Vr ± ev,, 


n 








pc 


10^ yr 


mas 


km s~^ 




km s~^ 


1 


2 


3 


4 


5 


6 


7 


8 


89825 


GL 710 


0.21 


1444 


51.12 ± 1.63 


-13.8 ±0.3 


4 


8.1 


70890 


Proxima Cen 


0.89 


27.4 


771.64 ± 2.60 


-25.1 ±0.9 






71683 


a Cen A/B * 


0.91 


28.4 


754.81 ±4.11 


-24.7 ±0.4 


1 




57511 


AC+79 3888 


1.0(3 


16.0 


186.86 ± 1.70 


-111.6 ±0.2 


5 


3.9 


87937 


Barnard's star 


1.15 


9.8 


548.31 ± 1.51 


-110.6 ±0.2 


5 


2.2 


27288 


GL 217.1 


1.27 


-862 


46.28 ±0.16 


24.7 ± 1.1 


2 


5.9 


54035 


Lalande 21185 


1.43 


20.5 


392.64 ±0.67 


-85.8 ±0.1 


7 


3.1 


26335 


GL 208 


1.54 


-500 


88.97 ± 1.02 


22.0 ±0.2 


5 


2.1 


57548 


Ross 128 


1.92 


72.4 


298.04 ± 2.30 


-31.0 ±0.2 


4 


1.1 


110893 


GL 860 A 


1.93 


90.0 


249.94 ± 1.87 


-34.0 ±0.1 


6 


2.3 


92403 


Ross 154 


1.96 


157 


336.72 ± 2.03 


-10.7 ±0.2 


5 


1.9 


38228 


HD 63433 


2.02 


1366 


45.45 ± 0.53 


-16.1 ±0.1 


6 


1.2 



Note, n is the number of catalogs used to calculate the mean V^; sy^ is the radial velocity error 
calculated from the external convergence (Gontcharov 2006); (*) the encounter parameters were 
calculated from component A. 



Table 3: Hipparcos stars from the solar neighborhood 30 pc in radius encountering with the solar 
orbit within d < 2 pc 



HIP 

1 


2 


SP 

3 


M/Mq 
1 


Rcf 

5 


lO''^ yr 
6 


7 


89825 


GL 710 


K7V 


0.4-0.6 


(1) 


0.311 ±0.167 


1447 ± 60 


70890 


Proxima Cen 


M5Vc 


0.15 ±0.02 


(2) 


0.890 ±0.019 


27.4 ±0.1 


71683 


a Ccn A/B * 


G2V/K1V 


2.0 ±0.1 


(2) 


0.910 ±0.012 


28.4 ±0.1 


57544 


AC+79 3888 


M4 


« 0.15 




1.059 ±0.013 


46.0 ±0.3 


87937 


Barnard's star 


sdM4 


0.144 


(3) 


1.148 ±0.006 


9.8 ±0.1 


27288 


GL 217.1 


A2Vann 


2.0 


(1) 


1.275 ±0.064 


-861 ± 40 


54035 


Lalande 21185 


M2V 


0.39 


(4) 


1.426 ±0.005 


20.5 ±0.1 


26335 


GL 208 


K7 


0.47 


(4) 


1.537 ±0.050 


-500 ± 8 


14754 


WD 0310-688 


DA 


0.5 




1.607 ±0.190 


-297 ± 29 


57548 


Ross 128 


M4.5V 


0.24 


(4) 


1.920 ±0.028 


72.8 ±0.7 


110893 


GL 860 A 


M2V 


0.56 


(4) 


1.929 ±0.028 


90.3 ±0.5 


92403 


Ross 154 


M3.5Ve 


0.23 


(4) 


1.959 ±0.030 


157 ± 1 


27887 


HD 40307 


K3V 




(5) 


1.974 ±0.035 


-413 ± 5 


38228 


HD 63433 


G5IV 




(5) 


2.038 ± 0.099 


1366 ± 20 



Note. The stellar mass estimates were taken from the following papers: (1) Garcia-Sanchez et al. 
(1999); Anosova et al. (1994); (3) Giampapa et al. (1996); (4) Mullari and Orlov (1996); and (5) 
Holmbcrg et al. (2007). If there is no reference, then a typical mass is given; (*) the sum of the 
masses is given. 
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Fig. 1. (a) Stellar trajectories relative to the Sun. The trajectories of the stars found here are 
highlighted by the thick lines: HIP 27887 (1), HIP 14754 (2), HIP 24186 (3), the pair of HIP 
104214 and HIP 104217 (4), HIP 99461 (5), and HIP 105766 (6). (b) The model trajectories 
of four stars computed by taking into account the random errors in the observational data (300 
realizations). The trajectories hatch the 3a confidence regions; the Oort cloud region is shaded. 




Fig. 2. Model trajectories of the star HIP 89825 (GL 710) relative to the Sun computed by taking 
into account the random errors in the observational data (300 realizations) . The trajectories hatch 
the 3cr confidence region; the Oort cloud region is shaded. 




Fig. 3. Model encounter parameters of the star HIP 89825 (GL 710) with the Sun. (a) The 
distribution of minimum distance dmin', (b) the histogram of encounter times tmin- 



